Spin Dephasing in the Extended Strong Collision Approximation 
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For Markovian dynamics of field fiuctuations we present here an extended strong collision ap- 
proximation, thereby putting our previous strong collision approach (Phys. Rev. Lett. 83 (1999) 
4215) into a systematic framework. Our new approach provides expressions for the free induction 
and spin echo magnetization decays that may be solved analytically or at least numerically. It is 
tested for the generic cases of dephasing due to an Anderson- Weiss process and due to restricted 
diffusion in a linear field gradient. 
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I. INTRODUCTION 



The understanding of spin dephasing is of paramount 
interest in all fields of nuclear magnetic resonance (NMR) 
sciences. In NMR spectroscopy it determines the line 
shape, in NMR imaging it is - besides longitudinal relax- 
ation - the major mechanism determining the contrast 
and contains morphological as well as functional infor- 
mation. 

The processes contributing to spin dephasing are re- 
lated to the spin environment. In biological tissues, for 
example, spin dephasing may result from dipole-dipole 
interaction of water proton spins with paramagnetic ions 
like Fe^^. Another cause is diffusion within inhomoge- 
neous magnetic fields generated by native or contrast 
agent induced susceptibility differences that are related 
to tissue composition and/or cellular and sub cellular 
compartments. In magnetic resonance imaging spin de- 
phasing in external gradient fields is exploited to get 
information about diffusion within biological systems. 
These diffusion sensitive imaging techniques are applied 
to study tissue anisotropy and restrictions of diffusion 
that are given by membranes of cells and sub cellular 
structures. 

Essential for dephasing of spins are the field fluctua- 
tions that induce the phase modulations. It is important 
to note that in biological tissues the relevant processes 
cover almost the whole range of time scales. For exam- 
ple, the dynamics of interactions of water proton spins 
with paramagnetic macromolecules as ferritin is so fast 
that it can be considered to be within the motional nar- 
rowing regime. On the other hand, dephasing of spins 
in magnetic field gradients around larger vessels is al- 
most coherent, i.e. it is in the static dephasing regime 
Hence, for biological applications it is important to 
obtain results from theory that are valid over the whole 
motion regime. However, in most cases this is not possi- 
ble analytically. 

Therefore, most efforts have focused on limiting cases. 
The motional narrowing limit is well investigated and 



a number of analytical results were obtained for it [0. 
The characteristic of this limit is that the mean phase 
shift induced by a field realization is much smaller than 
one, i.e. \Sf\ — r(Aw^)^/^ <^ 1, where the correla- 
tion time T gives the mean duration of some field re- 
alization, and (Acj^) is the variance of the inhomoge- 
neous field. The relaxation time is then obtained as 
I/T2 = T {IS.llP'). In the other limiting case, i.e. the 
static dephasing regime 

(r(Aw2)i/2 > 1)^ Yablonski |] 
derived analytical expressions for coherent dephasing of 
spins in inhomogeneous fields around magnetic centers 
like cylinders or spheres. Kiselev extended Yablon- 
ski's static dephasing approach by considering diffusion 
of spins within local linear gradients. However, this ap- 
proach requires that the diffusion length / during dephas- 
ing is within the linear approximation of the inhomoge- 
neous fields Lij{xq -\- I) k lj{xq) -\- dxio{xo) I. Note that 
expansion around the limiting cases by perturbation ap- 
proaches leads to divergences in the respective other lim- 
its. Therefore, the intermediate motion regime, i.e. al- 
most everything between the static dephasing and mo- 
tional narrowing limit, was in most situations accessible 
by simulations only Q. 

Recently, we used a strong collision (SC) approach 
to characterize spin dephasing in a particular situation: 
An inhomogeneous field around regularly arranged par- 
allel cylinders filled with a paramagnetic substance [Q, 
a model reflecting the capillary network of the cardiac 
muscle. The results agreed well with simulations ||^ over 
the whole dynamic range, and with experimental data 

§0- 

The basic idea behind the SC approach is to replace 
the original generator of the Markov process by a simpler 
one, the SC generator, that conserves particular features 
of the original process. In particular, by an appropriate 
choice of its parameter the SC process reproduces the 
correlation time of the field fluctuations induced by the 
original Markov process. There are several advantages 
of the SC approximation. First, it is is correct both in 
the motional narrowing and the static dephasing limits; 
thereby also the error in the intermediate regime is re- 
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duced considerably, when compared to perturbation ap- 
proaches. And second, it provides a simple expression 
for the magnetization decay which may be solved analyt- 
ically or at least numerically. 

However, the drawback of the SC approach was - up 
to now - that is not part of a systematic approximation 
to or an expansion of the original generator. Therefore, 
it was unclear how results could be improved beyond the 
SC approximation. The aim of this paper is to extend the 
SC approach and provide a framework for a systematic 
approximation. 

In the next section we will present a formal description 
of spin dephasing that will be the basis for our analyti- 
cal analysis. In Sect. Ill we will introduce the extended 
strong collision (ESC) approximation proper and show 
how it is used to describe free induction and spin echo 
decay. In Sect. IV we will apply it to two generic cases: 
spin dephasing induced by an Anderson- Weiss process ||] 
and by restricted diffusion in a linear field gradient. We 
will close the paper with a summary and a discussion of 
our results. 



II. FORMAL DESCRIPTION OF 
SPIN-DEPHASING 

We assume that dephasing of transversely polarized 
nuclear spins exposed to an external field is induced by 
randomly fluctuating magnetic perturbation fields with 
frequency Wi, where i is a discrete or continuous variable. 
The transition dynamics between two distinct states i 
and j is that of a stationary continuous time Markov pro- 
cess described by rates rji for the transition i j. The 
matrix R = (rji) as the generator of the Markov process 
conserves the probability to find a spin within one state, 
i.e. ru = — ^j^ifji- The eigenvalues I of R fulfill the 
condition I < where I = corresponds to the equilib- 
rium probability distribution. To simplify the notation 
we denote the normalized left and right eigenvectors of 
R as (/| and |/), respectively, with (/'|/) = Sin. 

The time evolution between t and t + dt of the trans- 
verse magnetization of spins in the state j (in polar no- 
tation rrij = rrij ^ — irrij y) results from the linear super- 
position of the transition and the precession dynamics, 
i.e. dtrrijit) — rji niiit) + i ujj mj{t). The precession 
within the external field was omitted since it only induces 
a constant offset of the frequency which may be gauged 
to zero. With the diagonal frequency matrix = (5ji Ui) 
one obtains for the magnetization \m) = {rrij) 

dt\m{t)) = (R + in) \m{t)) , (1) 

which is a generalization of the Bloch Torrey equation 
originally formulated for diffusing spins, i.e. R '--^ V^. 
In most cases it is reasonable to assume that the initial 
magnetization |to(0)) is proportional to the equilibrium 
probability distribution |0), e.g. when free diffusion is 
considered this would imply a homogeneous transverse 



magnetization. Equation (|l|) then provides the time evo- 
lution of the transverse magnetization (free induction de- 
cay) as 

\m{t)) = exp[(R -I- if2) t] |0) , (2) 

where the initial magnetization was normalized to one. 
The overall magnetization is then determined as 

M(t) = (0|m(t)) 

= (0| exp [(R + ir2) t] |0) . (3) 

The free induction decay as determined by Eqs. (||) and 
results from coherent and incoherent spin dephasing. 
The incoherent contribution is determined from spin echo 
experiments. In-plane polarized spins are rotated by a 
180° (TT-pulse) after a time t/2. This pulse transforms 
the original magnetization \m(t/2)) to its complex ad- 
joint \m*{t/2)) = exp[(R- iri) t/2] |0). This procedure 
cancels the coherent spin dephasing after the time t (echo 
time), i.e. the decay of magnetization at t is solely due 
to incoherent spin dephasing. The time course of the 
magnetization after the pulse, i.e. for times t' > t/2, is 
determined by 

\m{t')) = exp[(R + i n){t' - t/2)] exp[(R - iO)t/2]|0) , 

(4) 

i.e. the overall spin echo magnetization at the echo time 
t is 

Mssit) = (0| exp[(R + i n)t/2] exp[(R - iU) t/2] |0) 

^{m{t/2)\m*it/2)) . (5) 

Equation (^ relates the overall spin echo magnetization 
to the magnetization of the free induction decay. 

III. THE STRONG COLLISION 
APPROXIMATION AND ITS EXTENSION 

The analytical determination of the free induction de- 
cay according to Eqs. (^) is restricted to very few cases, 
e.g. free diffusion in a linear gradient or stochastic fluc- 
tuations between two precession frequencies. The idea 
of the strong collision approach and its extension is to 
replace the generator of the Markov process R by a more 
simple generator D that conserves specific features of the 
original dynamics. 

A. Strong Collision Approximation 

In many cases the stochastic fluctuations of the per- 
turbation fields occur on a much shorter time scale than 
spin-dephasing, i.e. the correlation time r of field fluc- 
tuations is much shorter than the relaxation time of the 
magnetization. For ergodic Markov processes one can 
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estimate that after a few r a spin has visited almost all 
relevant states with the equilibrium probability. On the 
other hand, there is only little change of the magnetiza- 
tion during this time interval. Therefore, spin dephas- 
ing in this situation can be described equivalently by a 
process in which the transition rate between two states 
i —f j is independent of the initial state. Consequently, 
the transition rate for i — > is proportional to the equilib- 
rium probability of the final state, po.j- Such a dynamics 
is referred to as strong collision dynamics. 
The generator D of this process has the form 



D = -A (id - Ho) 



(6) 



where IIo = |0)(0| is the projection operator onto the 
eigenspace generated by the the equilibrium eigenvector 
of R, and id is the identity operator. The factor A has 
to be determined self consistently. 

Since the starting point of the strong collision approx- 
imation is the observation that — in many cases of inter- 
est — the correlation of the stochastic field fluctuations 
appear on a shorter time scale than that of changes of the 
magnetization, only the long time behavior of the field 
fluctuations is of importance. This long time behavior 
is characterized by the correlation time of the two-point 
autocorrelation function of the field fluctuations (see also 
the appendix), 



C2(i) = {ioit)u{0)) 

= (0|J7exp(Rt)rj|0) 



(7) 



which is defined as 



T2 



dt 



C2W-C2(00) 
C2(0)-C2(00) 



(0|rj [exp(Rt) - Ho] f2|0) 
(0|O2|0) - (0|fl|0)2 



(8) 



Stochastic field fluctuations determined by the strong col- 
lision (SC) process should have the same correlation time 
as the original process, leading to the self-consistency 
condition 



ASC) 



(A) 



T2 



(9) 



It is easy to determine that the correlation time for the 
strong collision approximation is T2^'~^\\) — A^^, see the 
appendix Eq. (|C^), leading to 

A = r^-i . (10) 



B. Extended Strong Collision Approximation 

The extension of the strong collision approximation is 
based on a comparison with the spectral expansion of the 
original operator. 



R 



(11) 



where Zq = > Zi > . . . are the ordered eigenvalues of 
R and Ilj = is the projection operator onto the 

eigenspace corresponding to Ij . Since the time evolution 
operator is exp(R t) = X^jlo clear that the 

low order eigenvalues determine the long-time behavior, 
while higher orders dominate shorter and shorter time 
scales. A comparison with a rewriting of Eq. (^), 



D ^ /oHo - A (id - Ho) 



(12) 



(note that Zq — 0), shows that in the strong colli- 
sion approximation just the lowest order term of ( pT| ) 
is taken into account explicitly, while the contribution 
of the higher eigenvalues is approximated by the self- 
consistently determined parameter A. 

A natural extension, therefore, would be to take into 
account more low order eigenvalues explicitly, thereby in- 
creasing the accuracy of the description of the long-time 
behavior: 



Ij - A(id - n) 



(13) 



with n 



^^^Q IIj. A stochastic process generated by 
an operator D'„ in Eq. (^3|) will be referred to as sim- 
plified extended strong colHsion (ESC'„) approximation 
of order n. As before, the contribution of the higher 
eigenvalues is approximated by the parameter A, which 
is determined again self-consistently from condition (|^). 
Here it leads to 



c2(o)-e;=i 



C2(0)r2+E-=i^, 



(14) 



Note that 



with Wo, = (0|f2|i) and C2(0) = -^q^^. 
for n — > this equation becomes (|o|) again. 

However, there are several problems involved with an 
approximation based on ( p^ and ( p^ ) . Practically, an ex- 
act determination of the low order eigenvalues and eigen- 
vectors is possible only in special cases. Therefore, one 
has to deal with the problem that the eigenvalues and 
eigenvectors are known either only approximately or not 
at all. Furthermore, even with eigenvalues and eigenfunc- 
tions known, it turns out that the ESC approximation 
may be not applicable at all in certain situations: If the 
autocorrelation function of the field fluctuations is deter- 
mined fully by the eigenfunctions included in D', Eq. ( p^ ) 
is undetermined. In that case additional self-consistency 
requirements would be necessary for a better description 
of the intermediate time regime. 

Nevertheless, the above approach can be readily 
adapted to these situations. Equation ( p^ ) can be viewed 
as an optimized reduced description of the relaxation in 
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various subspaces of the original operator R. Such a op- 
timized description should also be possible for subspaces 
that are not eigenvectors of R. We can, therefore, set 



D„ = - ^ A, n, - A(id - n) 



(15) 



However, now the rates A,- 



1 , . . . , n are not eigen- 



values anymore, but have to be determined by additional 
self-consistency requirements, see below. Moreover, the 
Uj are not projectors onto the eigenspace of a particu- 
lar eigenvalue, but onto the spaces defined by arbitrarily 
chosen mutually orthogonal functions \ fj), j = 1, . . . ,n, 
with {fi\fj) = 5ij and (/j|0) — 0; i.e. the projectors have 
the form Uj = |/,)(/, | and H = Hq + Y."=i Natu- 
rally, one would try to choose the functions close to 
the eigenfunctions although it is not required for the 
extension to work. Another natural function space for 
example is based on polynomials in the frequency oper- 
ator r2, i.e. 



(16) 



where pi is some polynomial of degree i, the coefficients 
of which are chosen in such a way that the orthonormal 
relations are fulfilled. In the following we will refer to 
this base of functions as the fi-base. 

In analogy to Eq. ( [1^ ) a stochastic process generated 
by an operator D„ in Eq. ([l5| ) will be referred to as ex- 
tended strong collision (ESC„) approximation of order n. 
It is evident that the ESCq approximation refers to the 
strong collision approximation. 

We mentioned already that the rates Xj, j = 1, . . . , n, 
in Eq. (^5|) have to be determined now by additional self- 
consistency requirements. As it was with the SC ap- 
proximation, the aim of the ESC„ approximation is to 
approximate more closely the correlation of field fluctu- 
ations. This is achieved by considering also higher order 
correlation functions 

Cm {im—1 : ■ ■ ■ 1 ^1 ) — 



|r2exp(Rt™_i)r2...exp(Rti)rj|0) 



(17) 



Following the same arguments as for the strong collision 
approximation, the long time behavior of the Cm is of 
interest. In the same way as for the strong collision ap- 
proximation this should be characterized by some first 
order statistical moment which is obtained by integra- 
tion of the correlation function over tm-i, ■■■,ti. How- 
ever, direct usage of Cm is hampered by its non- vanishing 
asymptotic behavior: It is easily seen that from the rela- 
tion limf^^^oo exp(Rt,y) = Ho follows 



lim Cm{tm-1, ■■, il) — Cm-iy{tm-l, ■■, ii/+l) X 
— >oo 



which does not necessarily vanish. In the strong colli- 
sion approximation we avoided this problem by consid- 
ering the operator [exp(Rt) — IIo] instead of exp(Ri) in 
Eq. (^, i.e. we considered only the relaxational part 
of the stochastic process. When we perform the same 
replacement in Eq. ( |l7|) we obtain modified correlation 
functions c,„(t,„_i, . . . , ii) that we will call quasi cumu- 
lants (s. appendix). They vanish asymptotically for all 
t^. We now require that the generalized correlation times 
derived from these quasi cumulants. 



f H™ "'^ dt ' 

JO * Cm(0, . . 



,0) 



(19) 



are equal for the exact process and for the extended 
strong collision description. The relaxation rates are, 
therefore, determined by 

^(ESC)^-^^ Ai, . . . , A„) = Tm , 

TO = 2,4...,2n + 2, (20) 

which replace the single self-consistency condition (||). 
Note that in many systems the correlation functions 
Cmitm-i, ■ ■ ■ ,ti) vanish for odd values of m due to sym- 
metry. Therefore, we require the equivalence of relax- 
ation times in Eq. ( pO[ ) for even values of m only. Oth- 
erwise one has to determine the correlation times of the 
first n + 1 non-vanishing correlation functions. 

It is important to emphasize some properties of the 
ESC approximation. First of all, it usually does not re- 
duce to the ESC approximation when eigenfunctions are 
used for the projection operator; i.e. the Ai, . . . , A„ do 
not take on the numerical values of the corresponding 
eigenvalues, although they usually do approximate them. 
In the light of the problems with the ESC approxima- 
tions mentioned above, it will turn out that this is an 
advantage: the self-consistent determination of the relax- 
ation parameters A, and Ai, . . . , A„ according to Eq. ( [20| ) 
is more balanced than when only Eq. (^ is used, and 
gives rise to an improved approximation. Moreover, 
the self-consistency conditions ( pO| ) imply that both pro- 
cesses, the ESC process and the original Markov process 
have the same motional narrowing expansion of the t rans - 
verse relaxation, as it is shown in the appendix, Eq. 



(18) 



C. Transverse Spin Relaxation in the Extended 
Strong Collision Approximation 



In this section we will exploit the simple structure of 
the generator D„ to determine the time course of mag- 
netization. We will consider both: the free induction 
decay, i.e. the superposition of coherent and incoherent 
spin dephasing, and the spin echo decay, i.e. pure inco- 
herent spin dephasing. 
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1. Free Induction Decay 



2. Spin-Echo Decay 



In the extended strong collision approximation the gen- 
erator R of the free induction decay in the generalized 
Bloch Torrey Equation (|^) is replaced by the genera- 
tor D„ of Eq. (|lq). Instead of solving the propagator 
U(i) = exp((D„ -|-ir2)t) it is more convenient to solve its 
Laplace transform U(s) = (s — D„ — if2)~^ which may 
be expanded as 



U(s) = Uo(s + A) + Uo(s + A) A U(s) 



(21) 



where Uo(s) = {s — ifJ)^^ is the Laplace transform in 
the static dephasing limit (D„ = 0), and the operator A 
is defined as 



3=0 



(A-A,)m 



(22) 



where we set Aq = 0. We will now confine the opera- 
tors in Eq. ( pi] ) onto the subspace defined by the pro- 
jection operator 11 = X]J=o-'^i' Using the abbreviation 
O'^ := n O n for denoting any operator O confined to 
that subspace, we obtain 

U" (5) - U" (5 + A) + tjn {s + A) A^ tjn (5) , (23) 

where we exploited the fact that A = 11 A 11 and 
the idempotency of projection operators, i.e. 11 = 
n^. Equation (^) is of fundamental importance. It 
demonstrates that the relation (21) between the ESC- 
approximation and the static dephasing is also valid in 
the subspace [ |0), .., |/„) ]. This simplifies determi- 
nation of spin relaxation considerably since one only has 
to determine the (n-l-1) x (n-l-1) matrices ^ of the static 
dephasing limit and A, i.e. 

\5^{s) = (n - !!?(. + A) A^)-^ U?(,s + A) . (24) 

The Laplace transform of the overall magnetization de- 
cay M[„](s) in the extended strong collision approxima- 
tion has the form 



Mn(s) 



t^{s) |0) 



(25) 



For the special case of the strong collision approximation, 
ESCo, Eqs. (PI) and ^ result in 



Msdjs + A) 
1 - Msd[s + A) A 



(26) 



with Msd{s) = (0|Uo(s -I- A)|0) as the Laplace trans- 
form of the overall magnetization in the static dephasing 
regime. The time evolution M(t) can be obtained from 
Eqs (|2j) and (|2^) either by numerical inverse Laplace 
transform or using the generalized moment approach 
I0[ which allows a multi-exponential approximation. 



The spin-echo decay is obtained by inserting of the 
generator D„ into Eq. (||), i.e. 



dtM. 



SEAn 



it) 



(27) 



-XMsE.ln]{t) 
MsE ,[n]it) 
-AMsE,[„](t) 



(m(t/2)|A|m*(t/2)) 

(0|U(t/2)AU*(i/2)|0) 

(0|Un(V2)AU*n(t/2)|0) 



i.e. the spin-echo decay is expressed as a function of the 
spin-echo amplitude Mse, and the projection of the free 
induction decay onto the subspace [ |0), |/i), |/„) ], i.e. 
U'^(^)|0). This projection of the free induction decay is 
obtained from Eq. (24)by inverse Laplace transform, i.e. 
Un(t)|0) = /:-i(Un(s)|0)). The solution of Eq. is 



MsE,[n](t) = e 



-At 



rt/2 

1 + 2 d^e2^«(0|Un(^)AU*n(e)|0) 
Jo 

(28) 



1/T* = -\Ti{M{t))/t 

l/T2 = -HMsE)/t. 



M-i(M) 



dt M{t) 



3. Time Constants of Transverse Relaxation 



The free induction and the spin echo decay are usually 
described by the time constants Tg* and T2. However, 
there is no unique definition of these parameters. One 
definition of the relaxation times is 



(29) 



For the ESC decay one has to replace M by M[„] and 
Mse by MgE,[n]- This definition implies a dependence 
of relaxation times on t, except for single exponential 
decay. 

Another definition of relaxation times is based on the 
assumption that these constants provide the best single 
exponential approximation of magnetization decays, i.e. 
M{t) « e~*/^2 , Mse ~ e~*/^^ According to the mean 
relaxation time approximation the relaxation times are 
then the first long time moments of the decays |]l0|, i.e. 



^in case of degeneracy of the eigenvalues the matrix dimen- 
sion is the sum of the dimensions of the eigenspaces plus one 
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T2 ■.= fi-i{MsE) = / dtMsEit) 



(30) 



For a single exponential the mean relaxation time defini- 
tion and the definitions ( p9| ) give the same results. Ac- 
cording to definition ( ^0[ ) the relaxation times of the ESC 
decays can be related to their Laplace transforms as 



r; = M[„](o) 

T2 -M5E,[„](0) 



(31) 



The term M[„] (0) , which provides T2 , is obtained from 
Eq. (^5|). Applying some rules of Laplace transforms the 
term Mg^; r„i(0) giving T2 is obtained from Eq. (^ ) as 



This implies that the spin echo relaxation time is almost 
identical with the relaxation time of the absolute squared 
magnetization of the free induction decay, or vice versa 
that the free induction decay is almost irreversible. 

The dissipation-fluctuation-coherence relation ( |3^ ) 
takes a very simple form in the strong collision approxi- 
mation, when we assume that the overall magnetization 
decay is well approximated by a single exponential, i.e. 
M[o{{t) « e-*/^2-_ gj^(,g X = T^^, see Eq. (0), Eq. ^) 
reads 



To 



T2 
T2 



260 



(37) 



T2 =X-'+2 



E 

i=0 



(1 - A,/A) e. 



(32) 



where 



60 = 



^ / dz(0|Un(z)|0)(0|U*n(_z)|0) 



(33) 



is the mean relaxation time of the absolute squared 
overall free induction magnetization |M[„](t)p = 
|(0|U(t)n|0)|2, and for i > 1 



e,= / (0|Un(t)|/,)(/,|U*n|0)(t)) 
/o 

1 



27ri 



dz(0|Un(z)|/.)(/.|U*n(-z)|0) (34) 



are transit times describing the transient occurrence of 
the non-equilibrium components of the free induction de- 
cay U(t)|/i). Equation (|3|) relates T2 which describes 
the incoherent, i.e. irreversible, component of spin de- 
phasing to the stochastic field dynamics (A, Xi) and time 
constants of the free induction decay (6i), i.e. Eq. (^) 
is a dissipation-fluctuation-coherence relation. Note that 
the Eqs. ( p3[ ) to ( |3^ ) directly relate the time constants 
Qi to the Laplace transform of the free induction decay 
U'^(s) given by the fundamental Equation ([2^). 

^From Eq. ( p2| ) one can derive asymptotic relations 
for very fast and slow stochastic field fluctuations. Let 
e be some scaling parameter of D„, i.e. A, A^ e, then 
Eq. ( ^ ) reads in the static dephasing limit (e 0) 



To « A- 



(35) 



where we exploited that 8i(e) approaches its finite static 
dephasing limit. For very fast fiuctuations i.e. in the 
motional narrowing limit (e — s- cx)) one exploits that 
(/.|U(s)|0>/(0|U(.)|0) ^ 
demonstrates, i.e. one obtains 



as a power expansion 



T, « 2 Gn . 



(36) 



From Eqs. (37) follows that in the motional narrow- 



ing limit T2 ~ T2 holds whereas in the static dephasing 
limit of the strong collision approximation the relation 
T2 w T2 holds. 



IV. APPLICATIONS 



A. Anderson- Weiss Model 



The Anderson- Weiss model [g] is one of the rare ap- 
proaches - besides the ESC approximation - which de- 
scribes spin dephasing over the whole dynamic range of 
stochastic field fluctuations. The approach is suitable, 
for example, when dephasing is induced by spin inter- 
action with a great number of independently fluctuating 
perturbation fields in the spin environment. Then ana- 
lytical results are obtained for the free induction and the 
spin echo magnetization decay as 



M{t) = exp 
MsE{t) = exp 



f\t~0 C2{0d^ 
Jo 



t/2 



{t/2-0 C2i0d^ 



I {t~0 C2(C)rf? 

Jo 



(38) 



(39) 



where C2 is the modified two point correlation function 
(s. appendix). In this section we will first characterize 
the class of Markovian processes which fulfills the con- 
ditions of the Anderson- Weiss model. This leads to a 
generalized Bloch-Torrey equation according to Eq. (|l|) 
which is solved. Finally we compare the Anderson- Weiss 
model with its ESCq and ESCi approximation. 



6 



1. Markovian and Anderson-Weiss Dynamics 

The Anderson- Weiss approach is based on a Gaussian 
distribution of perturbation field frequencies lu. Even 
more important is the additional assumption that the 
stochastic phase accumulation of a spin cj) = w(^) 
also exhibits a Gaussian distribution. This latter condi- 
tion implies that the Greens function G{ujj,uJi, t), i.e. the 
probability that a spin initially precessing with frequency 
uji precesses at t with Uj, is also a Gaussian function in 
ujj^LOi with the condition 0(0;^, w^, 0) = 5{ujj — Ui). This 
implies that only nearest neighbor transitions rates are 
non- vanishing. Markovian processes in a continuous vari- 
able Lo with this property are described equivalently by a 
Fokker Planck Equation i.e. the probability density 
p{uj) satisfies 



dtp{uj,t) = 'Rp{u,t) 



F{uo))p{u,t) 



(40) 



where D{uj) is a ~ possibly tj-dependent - diffusion co- 
efficient and F{u)) is some driving force. Since the equi- 
librium probability density is a Gaussian function one 
obtains F(uj) = —c ■ uj with c > 0. The generalized 
Bloch-Torrey equation ([^) which determines the dynam- 
ics of magnetization as a superposition of precession and 
stochastic transitions then reads 

dtm(uj,t) = [d^D{uj){d^ + cu) +iuj]m{Lo,t) . (41) 



The derivation of the Eqs. (|4C|), ( f4l| ) is of fundamental 
importance since it states that a Markovian dynamics of a 
variable uj which satisfies the Anderson- Weiss conditions 
is equivalent to a diffusion process in this variable within 
a harmonic potential c uj"^ /2 and vice versa. Transfor- 
mation of variables uj — s- c^/^w and t — s- c^^^^t simplifies 
Eq. (|l|) to 

dtm{ij,t)^ [d^Piu;] [d^ + uj) + \ uj]m{uj ,t) (42) 

where we continue to denote also the transformed vari- 
ables as UJ and t and (3 — c?!"^ D is the transformed 
diffusion coefficient. In the following we will restrict con- 
sideration to the case of a constant diffusion coefficient f3. 
The left and right sided eigenfunctions of the transition 
operator R = pduj{duj+uj) arc the Hcrmite functions, i.e. 



with eigenvalues 



exp(-wV2)i/„(w) 
Hniuj) 



l^^ -n (i 



(43) 



(44) 



^From the definition of the Hermite functions and the 
operator intertwining relation [d^, (d^ -\- uj)] — 1 follow 
the recursive Equations 

\n + l) = -^d^\n) , |n - 1) = {d^ + u)\n) 



{n + l\ = {n\{d^ + uj) , {n - l\ = {n\d^ (45) 

n 

which also provide the normalization of eigenfunctions. 
The advantage of the Markovian formulation of the 
Anderson- Weiss model is that it does not only provide 
global parameters but also local ones, e.g. the time course 
of the magnetization with frequency uj. Straightforward 
application of the Eqs. ( ^5|) and some operator algebra 
provides the solution of Eq. (^2|) as 

m{uj,t) ^ exp[/39„(9^ +uj) + i uj]\0) 

= exp[-/3-H-f 

(27r)-i/2 exp[-l/2(co ~ i f]-\e>^' - 1))] (46) 

Integration over uj just gives the free induction decay of 
the overall magnetization 



M{t) = exp[-/3"4 + /3"2(1 - e"^^)] 



(47) 



which is just equivalent to the result of Eq. (38), since the 
2-point correlation function of Eq. (|4^) is C2 (t) — e^^* (s. 



appendix Eq. (El)). Insertion of this 2-point correlation 
function into Eq. ( |39| ) provides the spin echo decay as 



MsE{t) = M(i/2)2 exp[/?- 



-(3t/2 



If 



(48) 



Relaxation times of the free induction and spin echo de- 
cay were determined according to Eqs. (|30[). 



2. ESC approximation 

The ESC propagator is determined from the propa- 
gator in the static dephasing limit Uq = exp(iS7i) = 
[exp(iwi)], and the A-matrix of Eq. ( p2| ) both restricted 
either to the function space [ |0) ] for the ESCq or 
[ |0), 1/1) ] for the ESCi approximation Eq. (||). The 
special structure of the transition rate operator of the 
Anderson Weiss model implies that the base of eigenfunc- 
tions, Eqs. (^3|), is identical with the Sl-base, Eq. ([l^), 
i.e. p„(r2)|0) = |n). Hence, we will set in the following 
|/i)=Pi(f2) = |l) 

ESCq approximation The matrix element of the 
Laplace transformed propagator in the static dephasing 
limit required for the ESCq approximation is 



|Uo(s)|0) = y^e''/^ ceii{s/V2) 



(49) 



where cerf(z) = 1 — erf(z) is the complementary error 
function. The coefficient A which guarantees the self 
consistency condition Eq. (^ is determined form the 
Eqs. (C2,E2) (s. appendix) as 



(50) 
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ESCi approximation The matrix elements of the 
Laplace transformed propagator in the static dephasing 
limit required for the ESCi approximation in the Jl-base 
are that of Eq. (Egl) and 



(0|Uo(s)|l) = A^/:[(0|exp(int)O|0)] 

= (-i)/:[at(0|exp(im)|0)] 
= 7Vi(l-s(0|C/o(s)|0)) , 



(51) 



where the factor N generally is some normalization fac- 
tor with iV^ = (0|r2^|0), i.e. in the case of the Anderson 
Weiss model it is simply = 1. Consequently, using 
some elementary rules of Laplace transforms, one derives 
the other matrix elements as 



(52) 



(l|Uo(s)|0) = (0|Uo(s)|l) 
(l|Uo(s)|l)=7V2s(l-s(0|Uo(s)|0)) 



It has to be stressed that the Equations (|5l|) to (^ are 
generally valid for all ESCi approximations in the fJ- 
base. 

The coefficients A, Ai guaranteeing the sel f con sistency 
condition Eq. ( pO| ) are obtained from Eqs. (D4]E^) and 



Ai =/3 
A = 2/3 



(53) 



Relaxation in the ESCq and ESCi approximation 
The matrix tJj^(s) and the coefficients A and Ai deter- 
mine the Laplace transformed ESC propagator U'^(s) in 
Eq. (p^), which itself is the base for all other calculations. 
It directly provides Tj* when defined as the first moment, 
Eq. (|3l|), of the free induction decay Eq. (p5|). Insertion 
of Un(s) into Eqs. (||,||) provides according to Eq. (||) 
the relaxation time of the spin echo decay when defined 
as its first long time moment Eq. (|l|). Inverse Laplace 
transformation of U'^(s) gives the ESC propagator TJ^{t) 
which itself allows determination of the spin echo decay 

The relaxation time T2 := /i_i(M) of the Anderson- 
Weiss process is well approximated by the ESCq and 
ESCi approximation over the whole dynamic range of 
stochastic field fluctuations (Fig. 1). In the static dephas- 
ing regime all curves approach lim^^o mIi = ■\/2/7r. The 
successive approximation of the spin echo relaxation by 
the ESC approximation is seen from the magnetization 
decay curves (Fig. 2) and the curves showing the depen- 
dence of r2 obtained by either definition (Eqs. (|29|)-(|3^)) 
on the diffusion coefficient /3 as Figure 3 demonstrates. 
The latter curves all run parallel in the motional nar- 
rowing regime 

(ti((0|112|o))1/2 = p-i ^ 1) and exhibit 
a similar location of the maximum relaxation rate. To- 
wards the static dephasing regime (/3 —> 0) the rate of 
the Anderson- Weiss process declines less than the rates 
of the ESC processes. 



B. Spin Dephasing by Restricted Diffusion in a 
Linear Gradient Field 



1. The Exact Process 



Whereas dephasing of free diffusing spins in a linear 
gradient field can be treated analytically, only numeri- 
cal solutions exists for the restricted diffusion case 
On the one hand restricted diffusion in a linear gradient 
field provides a simple model to study principle features 
of spin dephasing by diffusion. On the other hand treat- 
ment of this problem is not only of academic interest as 
already mentioned in the Introduction. We will approxi- 
mate the free induction and spin echo decay of the global 
magnetization for the case of restricted diffusion by the 
strong collision approximation (ESCq) and its first exten- 
sion (ESCi). The ESCi approximation will be performed 
for both, in the f2-polynomial base, i.e. |/i) ~ JljO) and 
in the eigenfunction base, i.e. |/i) = 

We assume that the spins diffuse within an interval 
of size i in a linear gradient field uj(x) — g x. Re- 
flecting boundary conditions at x — ±L/2 imply that 
dxm{zLL/2,t) = 0. With D as the diffusion coefficient 
and R = E)[dl] the Bloch-Torrey Eq. (|l|) has the form 
dtm{x,t) = {D[dl] + igx)m{x,t), where the brackets [ ] 
denote that application of the operator is restricted to 
functions which fulfill the reflecting boundary conditions. 
Transformation of variables x ^ x/ L and t t g L re- 
sults in 



dtm{x, t) = (/3 \dl] + i x) m{x, t) , 



(54) 



and vanishing derivatives at the edges of the unit interval 



a^m(±l/2,t) EE 



(55) 



with the diffusion coefficient /? = D/{gL^). We continue 
to denote also the transformed variables as x and t to 
reduce the number of symbols. When the initial mag- 
netization is proportional to the equilibrium probability, 
i.e. m{x,Q) = 1, the Laplace transform m{x,s) of the 
local magnetization decay satisfies 



(P [d^] + i x)m{x, s) 



-1 



(56) 



The Equation ( ^^ was solved numerically. Integration of 
the result over the unit interval provided the free induc- 
tion decay of the overall magnetization, and application 
of Eq. (H) on the result gave the spin echo decay. When 
the spin echo relaxation time was defined as the first sta- 
tistical moment of the magnetization decay Eq. (|3^) was 
applied. For determination of Tg*, defined as the first 
moment of the free induction decay, Eq. ( |56| ) was solved 

numerically, and integration Xi(y2 ^•'^ "^(2^! ^) — M{s) 

gave T* = M(0) = n-i{M). 
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2. ESC approximation 

The determination of the ESCq and ESCi approxi- 
mation is completely analogous to that for the Ander- 
son Weiss model, except that the ri-polynomial and the 
eigenfunction base are not identical. 

ESCq approximation The equilibrium function for 
the restricted diffusion within the unit interval is 



(l|Uo(.)|l) 



|0> = 1 



i.e. one obtains 



|Uo(s)|0) =i In 



i/2 



■i/2 



(57) 



(58) 



The self consistency condition for the strong collision ap- 
proximation (|9[) determines the parameter A as (s. ap- 
pendix Eqs. (^2]) ( pTl| )) 



A = 10/3 



(59) 



Insertion of the results of Eqs. (||) and (|g) into Eq. (|2|) 
determines the Laplace transformed propagator in the 
ESCo approximation (0|U(s)|0) from which T2, T2 and 
spin-echo decay curves are obtained. 

ESCi approximation in the ^.-polynomial base: The 
lowest order function besides the equilibrium state in the 
fi-polynomial base has the form 



i/i> = mn 

= 2\flx 



(60) 



The matrix element (|5^) and the Equations (|51|) to 
(^) then directly provide the static dephasing operator 
Uo(s) = (s — i x)^^ in the [|0), |/i)] base. The parameters 
Ai,A of the ESCi approximation are determined from 
t he s elf consistency con di tion (po|), i.e. with Eqs. (D4), 
(pTl|)and and Eqs. (p7|),(pT7|) one obtains 



Ai = 10/3 
_ 443520 

^-^9or^ 

w 49.83 /3 



(61) 



Development in the eigenfunction space: The normal- 
ized non equilibrium eigenfunctions of the restricted dif- 
fusion operator are 



\v) = \pi sin(^7r x) for = 1, 3, . 
= \pi co^iy-K x) for = 2, 4, , 



(62) 



Since [9^] is a symmetric operator left and right sided 
eigenfunctions are identical. With |/i) = |1) and z = 
7r(l/2 -I- i s) one obtains 



(0|Uo(s)|l) = V2[sinh(^s)Ci(e) + 
i cosh(7rs)Si(0] 



-2 arctan(2s)— 
[i cosh(27rs)Ci(0+ 

sinh(27rs)Si(C)] 



(63) 



where Ci and Si denote the integral cosine and integral 
sinus function respectively. The parameters Ai A in the 
eigenfunction base are determined similarly as in the f2- 
base (s. appendix) and one obtains 



Ai 
A 



9.89/3 
41.6/3 



(64) 



Relaxation in the ESCq and ESCi approximation 
Figure 4 demonstrates the first long time moment of 
the free induction decay, that of the strong collision 
approximation, and its first extension for both bases 
as a function of the diffusion coefficient /3. All curves 
show the same asymptotic behavior in the static de- 
phasing (/3 — > 0) and in the motional narrowing limit 
(ri(0|a;2|0) = 1/120 < 1). Furthermore the bet- 
ter approximation by the ESCi curves compared to the 
ESCq curve in the intermediate motion regime is evident. 
There is no significant difference between the ESCi ap- 
proximation in the eigenfunction and in the H-polynomial 
base. 

The spin echo magnetization decay is shown in Fig. 5. 
Especially in the long time behavior near the static de- 
phasing regime, the ESCi curves either in the eigenfunc- 
tion space or in the Jl-space demonstrate a better ap- 
proximation than the ESCq curve. This is also reflected 
by the dependence of spin echo relaxation rate I/T2 on 
the diffusion coefficient (Figs. 6,7). When defined by the 
echo time (Eq. (p9|)) the ESCq and ESCi curves run par- 
allel with the curve obtained for restricted the diffusion 
dynamics for short echo times. For longer echo times and 
decreasing diffusion coefficients the ESCi curve provides 
a better approximation. Again as for the free induction 
decay there is no significant difference between ESCi ap- 
proximations in the eigenfunction and that in the fi-base. 



V. SUMMARY AND DISCUSSION 



Analytical results on transverse spin relaxation due 
to stochastic phase modulation exist mainly for limiting 
cases, like the motional narrowing and the static regime. 
Perturbation approaches are only valid close to their re- 
spective limits, and they diverge as one tries to extend 
them towards the opposite motion regime. Particularly 
the intermediate motion regime cannot be described re- 
liably by such a treatment. 

We choose a different approach. Our aim was to ap- 
proximate the dynamics, assumed to be Markovian, by 
a more simple one that conserves specific features of the 
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original. The starting point was the strong coUision ap- 
proximation that assumes the transition probabihty 
between two states being independent from the initial 
state, an approximation that holds when spin dephas- 
ing occurs on a time scale significantly longer than the 
stochastic phase modulations. Hence, all states perpen- 
dicular to the equilibrium state relax with the same ex- 
ponential factor that is determined self-consistently by 
comparison with the field fluctuations. 

Note that the motional narrowing limit as well as the 
static dephasing regime are described correctly by this 
approximation. Consequently, the error in the interme- 
diate motion regime is already less than it would be by 
perturbation approaches of a comparable low order. Nev- 
ertheless, there is still room for improvement. Also, one 
would like to have higher order approximations that can 
be used to check the quality of low order descriptions. 

A systematic extension of the strong collision ansatz is 
to include the relaxation of states of an appropriate fi- 
nite function base explicitly. We require that correlation 
times of original and approximate dynamics are identical 
to a certain order. This self consistency condition assures 
that both dynamics have the same motional narrowing 
expansion of spin dephasing. As it was already in the 
strong collision ansatz, spin dephasing is asymptotically 
identical for both dynamics in the limit of the static mo- 
tion regime. 

The finite function base of the ESC„ approximation 
may be given by the first n ordered eigenfunctions of the 
generator of the original phase modulations. Obviously, 
then the ESC generator directly reflects the dynamics of 
original generator up to a time scale corresponding to 
the n-th eigenvalue. For practical applications the ESC 
approach within an eigenfunction space may be a safe 
way to approximate spin dephasing. However, when the 
determination of the eigenfunctions is tedious, the appli- 
cation of the abase [|0), - rj|0), I/2) - fl^\0), . . .] 
may be more appropriate, at least for the ESCi approxi- 
mation. Within the fi-base the determination of the two- 
and four-point correlation times (s. appendix) and the 
propagator in the static motion regime, Eq. (|52|), is con- 
siderably simplified. 

The mechanism by which the ESCi approach in the Q- 
base works becomes evident by the following considera- 
tion: terms of the motional narrowing expansion Eq. (Bl) 
may be interpreted as repetitive interactions of the spin 
system with the inhomogeneous field $1 and intermedi- 
ate evolution with the free propagator exp(Rti). In the 
motional narrowing limit, one obtains from Eq. (B7) 



1/1)7 (/il/i) = 1- Equation ( p^ j implies that in the mo- 
tional narrowing limit the long time behavior of spin de- 
phasing solely depends on the free propagator related re- 
laxation of the state |/i), i.e. this state remains the only 
relevant one. Hence, it is obvious that in the interme- 
diate motion regime an ESCi approximation including 
the state |/i) in its generator is superior to the ESCq 
approximation. 

Within the function base the propagator of spin de- 
phasing is directly related to the propagator of spin de- 
phasing in the absence of stochastic phase modulations. 
This specific feature of the ESC dynamics tremendously 
facilitates the actual determination of spin dephasing for 
the following reasons: (i) in many cases the propagator 
in the static motion regime (which is an average phase 
factor) may be determined analytically or at least nu- 
merically; (a) the determination of the propagator from 
that in the static motion regime is self contained within 
the base, i.e. it is obtained from a combination of finite 
dimensional matrices. 

The two lowest order ESC approximations were ap- 
plied to two generic models: spin dephasing in the An- 
derson Weiss model, i.e. Gaussian frequency distribution 
and Gaussian transition dynamics, and dephasing by re- 
stricted diffusion in a linear frequency gradient. The 
reason for this choice was that - besides their generic 
character - these models allow either an analytical (An- 
derson Weiss) or, at least, a simple numerical treatment 
(linear gradient) of magnetization decay. These features 
are helpful to prove the ESC approach. For the An- 
derson Weiss model we determined the corresponding 
Markov generator of the phase modulations, which - to 
our knowledge - was done here for the first time. 

For both generic models the subsequent improvement 
by ESC„ approximations of dephasing parameters and 
magnetization decays could be demonstrated. One of our 
next aims will be the application of the ESC approach to 
more realistic scenarios. 

In closing, we would like to emphasize that the ESC 
approach is actually not limited to spin dephasing only. 
It can be applied, in principle, in any situation where the 
time behavior of complicated observables of stochastic 
processes is of interest. In each case, however, an appro- 
priate function base has to be chosen, corresponding to 
the r2-base for spin dephasing. 
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1/T2=C2(0) 



dt (0|nexp(Ri)J^|0) 

POG 

\n'\0) / (/i|exp(Rt)|/i) 



(65) 



where the factor (0|J1^|0) is due to the normalization of 
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APPENDIX A: AUTOCORRELATION 
FUNCTIONS AND QUASI-CUMULANTS 



APPENDIX B: MOTIONAL NARROWING 
EXPANSION 



The general n-point autocorrelation function of 
stochastically fluctuating fields cuj is defined as 



Cn{tr, 



,ti) 



n-l 



n-1 



E 



(Al) 



v=0 



where ^(cjj^ j , ^"^j^'^ t,; tjjj , ti; Wjo, 0) is the probabil- 
ity to find at t = the frequency Wj^, at i = ti the value 
Wjj,... and dXt = 'Y^^^i the frequency tjj^ ^.When the 
stochastic dynamics is determined by a Markov process, 
this probability can be factored into transition probabili- 
ties between sequential states i — > i -|- 1 after the interval 
ii+i and the initial {t = 0) probability distribution, i.e. 

T!-l 



i=l 

n-l 



(A2) 



i=l 



The transition probabilities after the interval ti are the 
matrix elements of evolution operator exp(Rii). Since 
the dynamics is assumed to be stationary the initial prob- 
ability p{ojjo , 0) is the equilibrium state probability dis- 
tribution, i.e. we can rewrite Eq. (Al) 



C„ = (0| r2exp(m„_i)0....exp(mi)r2 |0) , (A3) 

where fl — (ujjdj^k) is the diagonal frequency matrix. 
A modification of the correlation functions occurs if one 
exchanges the evolution operator exp(Rt) with the oper- 
ator exp(Ri) — IIo, where IIo — |0)(0| is the projection 
operator onto the equilibrium state space. This modified 
evolution operator describes the relaxation of observables 
minus their equilibrium state values. The modified auto- 
correlation functions will be denoted as quasi cumulants 
and they are then defined as 



c„ = (0| f2[exp(Ri„_i)-no]f2.. 
....[exp(Rti) -no]f^ |0) 



(A4) 



The Laplace transform of the correlation function in 
Eq. (A4) has the form 



c„(s„-i, ..si) 



1 

51 -R 



s„-i 

^ n 



n.... 

n |0) . (A5) 



This Laplace transform allows the determination of tem- 
poral moments of the normalized autocorrelation func- 
tion Cn{tn-i, ..,ii)/c„(0, ..,0) as the generalized correla- 
tion times 



c„(0,..,0)/c„(0,..,0) 



(A6) 



The motional narrowing expansion is a perturbation 
approach to determine the overall magnetization M{t) - 
or its Laplace transform M{s) - in terms of powers of the 
fluctuating fields fl. It is based on the assumption, that 
the stochastic fluctuations are more rapid than the pre- 
cession frequencies of the perturbation fields (motional 
narrowing limit). We will present a general relation be- 
tween the relaxation of the magnetization and the corre- 
lation of the field fluctuations that contains the motional 
narrowing limit as a limit case. The first step is to ex- 
pand the Laplace transform of the overall magnetization, 
Eq. d), in n, i.e. 



Mis) 



1 



|0> 



' ' s - R - irj 

(0| (s - R)"^ + i(s - R)-ifl(s - R)-^ 
+i\s - R)-iO(s - R)-ir2(s - R)-i -f- 



•|0) 



-2 ;2 



i (0|r2|0) 



i^ (0|f2(s- R)-ifJ) 



-i^l + .-i f;i'^a(s,s,..,s)') , 



(Bl) 



where C^, are the Laplace transformed n-point corre- 
lation functions of Eq.(A3). To avoid singularities at 
s = it is better to consider M~^{s). When we set 
q = s^^ J2T=i '^^ ^v{s, s, .., s) one obtains 

(ly-l)x 

M"i(,s) = s (l + Y,{^l)P qp\ 

oo oo 

= s-^i'^a + s-^ i^'+'^aia. + ... 

na„ + ... 



(B2) 



(-1)'' ^ 

iyi,..,iJp = l 



oo p 



m—1 



Rearrangement of terms of equal order in SI provides 

oo 

M-\s) = s-Y'^^ K, , (B3) 

i=i 

where the coefficients Kj have the form 
Kj = Cj s ^ ^ Ciy^ -t- . . 
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ui+..+Vf,=j m=l 



(B4) 



A comparison of t his s um w ith the modified correlation 
functions Cj, Eqs. (A4) and (A5), shows that 



(B5) 



i.e. one obtains 



(B6) 



0-1) X 



Equation ( [B6D expands the relaxation of the magneti- 
zation in terms of correlation functions to an arbitrary 
order. The long time behavior of is determined by 
the Laplace transform in the range of small s, i.e. in this 
range the relation 



c(0,0..,0) 

Cj(0,0..,0) 



(B7) 



is valid. The series in Eq. (B7) contains terms of mag- 
nitude < (0|r2''|0)/Z'^, where I denotes non vanishing 
eigenvalues of R. The latter determine the fluctuation 
frequency. In the motional narrowing limit these fluctu- 
ations are much higher than the precessing frequencies 
(0|r2''|0)/Z'^ <C 1, i.e. after normalization of the average 
perturbation field (0|r2|0) to zero, i.e. c\ = 0, M{t) is 
given a single exponential decay with the well known re- 
sult for the transverse relaxation rate as ^ = ^2(0) = 
n (0|f22|0) . 



APPENDIX C: QUASI CUMULANTS IN THE 
STRONG COLLISION APPROXIMATION 



In the strong collision (ESCo) approximation the quasi 
cumulants take a very simple form. Insertion of the gen- 
erator D = — A(id — Ho) into Eq. (A5) results in 



n — 1 

= n — -T cn^i(O'-o) , (CI) 

i.e. the quasi cumulant is a product of single exponen- 
tial f unc tions e~^*' and the generalized correlation times, 
Eq. ([A^), are all identical namely 



(C2) 



APPENDIX D: QUASI CUMULANTS IN THE 
EXTENDED STRONG COLLISION 
APPROXIMATION 

We restrict ourselves here to the ESCi Approximation 
and determine the correlation functions and generalized 
relaxation times for the fJ-polynomial base only in order 
to show the principle. Extensions to higher order approx- 
imations and to other function bases are straightforward, 
although they may be more tedious to calculate. 

For ESCi the generator of the stochastic field fluctua- 
tions has the form D — — Ai Hi — A (id — Ho — Hi). We 
will determine only the Laplace transforms of the 2- and 
4-point correlation functions since the 3-point correlation 
functions vanishes in the models we consider. According 
to Eq. ( A5) the determination of the correlation functions 
requires of the operator 



1 



s - D s 



1 



s + Ai 



s + A 



(id -Ho -Hi) 



(Dl) 



Assuming that the average frequency vanishes, i.e. 
(n) = (0|r2|0) = 0, which can always be achieved 
by normalization., ^ r2|0). Hence, the projector 
III = 1/1) (/i I takes the form 



Hi 



(D2) 



For the Laplace transformed two point correlation func- 
tion one obtains then 



.{ESCi, 



(s)-(Ai+.)- 



Ifl^lO) , 



(D3) 



i.e. the two point correlation function exhibits a single 
exponential decay with relaxation rate 



(BSCi) _ ,-1 
T2 — Ai . 

The 4-point correlation function is 

C\ (S3,S2,Si) = 

1 



1 



1 



(s3 + Ai)(si -f Ai) \S2 + Ai S2 + A 
1 



(D4) 



)|f22|0) 



(S3 + Ai)(s2 + A)(si -FAi) 
(D5) 

The four point correlation time is then determined as 

1/3 



SESCi 



1 / 1 



1 



|J^3|0)2 



Af VAi a ) (0|n2|o)((0|n4|o) - (0|n2|0)2) 



1/3 



-I- 
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(D6) 



2. Restricted Diffusion in a Linear Gradient Field 



In the case of the Anderson Weiss model and for the re- 
stricted diffu sion Unear gradient one has (0|11'^|0) = 0, 
i.e. Eq. (D6) simphfies to 



(ESCi) _ 3 



1 



(D7) 



APPENDIX E: RELAXATION TIMES OF 
CORRELATION FUNCTIONS IN THE MODELS 



In this final appendix we will determine the general- 
ized relaxation times of stochastic field fluctuations up 
to the 4-th order for the generic models we discuss in the 
main text. 



1. Diffusion in a Harmonic Potential 



When field fluctuations result from diffusion in a har- 
monic potential according to Eq. ( ^0| ) and spin dephas- 
ing is described by Eq. (^), the corresponding Laplace 
transformed 2-point correlation function, Eq. (A5), is 



C2(s) 



1 



' UJ- 



1 



^|0) 



+ /3 ' 



(El) 



where we applied the operator properties of d^j-, (du + uj) 
according to Eqs. p5|). This result shows that the two 
point correlation function exhibits a single exponential 
decay. Since C2(0) ~ (0|a;^|0) — 1 one obtains 



T2 = 



(E2) 



Since C3(t) vanishes, the next relevant correlation func- 
tion is c^it). Similarly, one obtains for its Laplace trans- 
form 



C4(S3, S2, Si) 



is3 + P)is2+2p){si+l3) 



(E3) 



And since 04(0) = (0|w*|0) — 2 the corresponding corre- 
lation time is 



T4 



(E4) 



In this s ecti on we will determine the correlation times 
r„ of Eq. ( A6) for n = 2, 4 for 1-d restricted diffusion of 



spins in a unit box in which they are affected by a linear 
gradient field. In dimensionless parameters one obtains 
for the generator R — f] [9^] where f3 is the dimensionless 
diffusion coefficient. The brackets denote that the appli- 
cation of the operator 9^ is restricted to functions with 
vanishing derivative at x — ±1/2 (reflecting boundary 
conditions). The frequency operator is Jl = a; and the 
equilibrium state eigenfunction is |0) = 1. 

The determination of correlation times r„ requires the 
calculation of t he L aplace transformed quasi cumulants 
c„(0, ..,0), Eq. (|A6|). The definition of this quasi cumu- 
lants according to Eq. (A5) shows the need of recurrent 
determination of terms of the form 



lim[-(s-R)-i+,s-ino]5, 
with some function g, i.e. / fulfills 
R/ - (1 - Ho) g , 



(E5) 



(E6) 



i.e. with R = [d^] Eq. ( [Eq ) becomes a second order dif- 
ferential equation. The application of the inverse second 
order differential operator leaves in general 2 integration 
constants. One may be determined from the reflective 
symmetric boundary conditions; however, a further con- 
dition is needed to get the second constant. A spectral 
decomposition R in Eq. (E5) shows HqJ = 0, i.e. we 
obtain as a further condition 



1/2 



dx f{x) = . 

-1/2 

Two point correlation time: We define 



(E7) 



/i(x) = lim(-(s-R)- 

s— ►0 



s-ino)x |o) 



= {x^/&-x/%) 



(E8) 



which obviously fulfills the reflecting boundary condition 
at a; = ±1/2 and Eq. (E7). Hence, the two point Laplace 
transformed quasi cumulant is 



C2(0) 



x lim(-(s - R)"i + s-^Uo) x |0) 

soO 



1/2 ^ ^ 

dx xfi{x) — (3^^ 

-1/2 



120 



and with 

C2(0) = (0|x2|0) 
one obtains 

T2 = 



1/2 



dx X 



1/12 



-1/2 



(E9) 



(ElO) 



(Ell) 
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Fo ur point correlation time: Iterative application of 
Eq. (E5) defines 



/2(x) = lim(-(s - R)-i + s-^Uo) xh{x) 



(E12) 



which, according to Eq. ([Eq), fulfills 

[dl\f2{x) = {l-U^)xh{x) 



xfi{x) 



1/2 



1/2 



dx xfi (x) 



xfUx) + 



(E13) 



Insertion of /i, Eq. (|E8|), and considering the reflective 
boundary conditions and the condition (|E7|) yields 



180 96 240 161280 



37 \ 



(E14) 



Similarly to the procedure above we could determine a 
function /3(x) but instead we use a different approach 
which exploits the symmetry of eigenfunctions of the op- 
erator R = The Laplace transformed 4-point quasi 
cumulant is 



=/i(^) 



C4(0,0,0) 



(0| X lim 

1 



1 



1 

5^ + -no 

S3 - R S3 



lim 

S2^0 



X lim 



R + - Ho 

S2 - R S2 



1 



In. 



1/2 
-1/2 



Si — R si 
fi{x)xf2{x) 



79833600 ' 



and with 

C4(0,0,0) = (Olx-^IO) - (0|a;2|0)2 
we finally have 



1 

180 



7"4 



443520 
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FIG. 1. Relaxation time T2 (defined as the first long time 
moment ^-1) of the free induction decay in the Anderson 
Weiss model (AW) and its ESCo and ESCi approximation as 
a function of the diffusion coefficient /3. 



FIG. 2. Spin-echo magnetization decay in the Anderson 
Weiss model (AW) and its ESCo and ESCi approximations for 
three different diffusion coefficients /3. Note: for the diffusion 
coefficient close to the motional narrowing regime (/3 = lO'''^) 
the Anderson Weiss and the ESC curves almost run paral- 
lel. Therefore for the clearness of the figure, only the An- 
derson Weiss curve is shown. In the intermediate motion 
regime 13 = 1 the original and the ESCi curve still run parallel 
whereas the ESCo approximation already shows a moderate 
deviation. Towards the static dephasing regime (/3 = 10""'^) 
the successive improved approximation of the Anderson Weiss 
curve by the ESC curves is evident. 



FIG. 3. Dependence of the spin echo relaxation time T2 
on the diffusion coefficient P for the Anderson- Weiss model 
(AW) and its ESCo and ESCi approximation: (a) The re- 
laxation time was defined according to Eq. (E9) by the 
echo time t, and (b) as the first long time moment /^_i of 
spin-echo magnetization decay according to Eq. (|3o|). The 
Anderson Weiss curves and the corresponding approxima- 
tions converge as /3 approaches the motional narrowing regime 
ri((0|n^|0))^/^ = /3"^ < 1. When defined by the echo time 
(a) the T2 curves of the Anderson Weiss model and its ap- 
proximations all run parallel for the short echo time {t = 1). 
With increasing echo time {t = 3, 6) the successive ESC ap- 
proximation becomes evident. 
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FIG. 4. Relaxation time of the free induction decay T2 of 
spins diffusing within a linear field gradient in the unit inter- 
val as a function of the diffusion coefficient /?. T2 is defined 
as the first long time moment u„i and obtained from the 
Bloch Torrey (BT) Equation (|j). The ESC approximations 
are shown. The ESCi approximation was determined for the 
eigenfunction (ef) and the fi-base. 
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FIG. 5. Spin-echo magnetization decay for restricted dif- 
fusion within a linear field gradient in the unit interval as 
obtained from the Bloch Torrey (BT) Equation (|^. Three 
diffusion coefficients /3 are considered. The ESC approxima- 
tions in the different diffusion regimes are demonstrated. The 
ESCi approximation was obtained for the eigenfunction (ef) 
and the fi-base. 



FIG. 6. Spin echo relaxation time T2 as a function of the 
diffusion coefficient (5 for restricted diffusion in the unit inter- 
val and the corresponding ESC approximations. The labeling 
of the curves is as in Figure 4. The relaxation time was de- 
fined by the echo time t according to Eq. jiol ) . 

FIG. 7. Spin echo relaxation time T2 as a function of the 
diffusion coefficient /3 for restricted diffusion in the unit inter- 
val and the ESC approximations. Labeling is as in Figure 4. 
The relaxation time is defined as the first long time moment 
;x_i of the spin echo decay Eq. (^). 
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